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Abstract 

In Part I soliton solutions to the ABS list of multi-dimensionally consistent 
difference equations (except Q4) were derived using connection between the Q3 
equation and the NQC equations, and then by reductions. In that work central 
role was played by a Cauchy matrix. In this work we use a different approach, we 
derive the iV-soliton solutions following Hirota's direct and constructive method. 
This leads to Casoratians and bilinear difference equations. We give here details for 
the H-series of equations and for Ql; the results for Q3 have been given earlier. 

1 Introduction 

The analysis of integrability for discrete systems is now in active development. Discrete- 
ness introduces many complications in comparison with continuous integrability, mainly 
due to the lack of Leibniz rule for discrete derivatives. This also includes the fact that 
space-time itself has many discretizations. 

The basic philosophy and definitions have been given in Part I p], here we repeat only 
some essential ingredients. The underlying space is formed by a Cartesian square lattice 
and the dynamical equation is denned on an elementary square of this lattice. (There are 
other possible settings but even this case has not been fully analyzed.) As for the defi- 
nition of integrability we choose "Consistency-around-the cube" (CAC) which is further 
explained in Section 12.11 With this choice of integrability many nice properties follow, 
including the existence of a Lax pair. Furthermore, with mild additional assumption one 
can classify the integrable models [2] and the 'ABS list" is surprisingly short. 

In this paper we construct multi-soliton solutions for the H-series of models in the 
ABS-list, as well as for the Ql model. Our approach is constructive for each model and is 
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based on the approach of Hirota, as explained in Section 12.21 Some work on this problem 
already exists for the two highest members of the ABS-list, but the Q4 results in [3 J is far 
from explicit and the result for Q3 [4] relied on a specific association with the NQC-model, 
which has been further elaborated in Part I. 

We hope that the detailed analysis of the simpler models, for which everything can be 
made systematic and explicit, will provide further understanding on the soliton question 
and that it could be used for other models as well. 

In the next section we will give the general background for our approach and then in 
subsequent sections we will go through the models HI, H2, H3 and Ql. 



2 Generalities 



2.1 Multidimensional consistency and the ABS list 

As in Part I we only consider quadrilateral lattice equations defined by a multi-linear 
relation on the values at the four corners of an elementary square of Cartesian lattice, see 
Figure l2TTl (a). Given a base point u nm = u we indicate the shifts in the n,m directions 
by a tilde or a hat as described in Figure 12. 1[ 

Multidimensional consistency is an essential ingredient in the construction of soliton 
solutions, the new dimensions standing for parameters of the solitons. The consistency 
is defined as follows: We adjoin a third direction and construct an elementary cube as in 
Figure l2Tl7 b). the shifts in the new third direction are denoted by a bar, e.g., w n ,m,ifc+i = u. 
On the base square we have the equation Q(u,u,u,u;p,q) = 0, and we now introduce 
the same equation with different variables on the sides and on the top, thus we will have 
altogether the following set of equations 

Q(u, u, u, u;p, q) = 0, Q(u, u, u, u;p, q) — 0, (2.1a) 

Q(u,u,u,u]p,r) = 0, Q(u, u, u, u;p, r) = 0, (2-lb) 

Q(u,u,u,u;q,r) = 0, Q{u, u, u, u; q, r) = 0. (2.1c) 

Now considering Figure l2~7fl (b) we can have initial values given at black circles (u, u, u, u), 
and in order to compute the remaining four values we have six equations. We can use the 




Figure 1: (a): The points on which the map is defined, and (b): the consistency cube. 
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LHS equations in ( 12 .lap . ( I2.1bl) . (12. left to compute u, u, u, respectively, and this leaves 

the three RHS equations from which we should get the same value for u, and this implies 
two consistency conditions. 

In [2] a classification of quadrilateral lattice equations was performed following the 
above definition of integrability, with two additional requirements on the equations: sym- 
metry and the so called "tetrahedron property", which means that the computed value 

for u should only depend on u, u, u and not on it. The ABS list is the following: 

(u — u) (u — u) + q — p = 0, (HI) 

(u — u) (u — u) + (q — p) (u + u + u + u) + q 2 — p 2 = 0, (H2) 

p(uu + uu) — q(uu + uu) + 5(p 2 — q 2 ) = 0, (H3) 

p(u + u)(u + u) — q(u + u)(u + u) — 5 2 pq(p — q) = 0, (Al) 

(q 2 — p 2 )(uuuu + 1) + q{p 2 — l){uu + uu) — p(q 2 — \){uu + uu), (A2) 

p{u — u) (u — u) — q(u — u)(u — u) + 5 2 pq(p — q) =0, (Ql) 



p{u — u){u — u) — q{u — u){u — u) 

+pq{p — q){u + u + u + u) — pq{p — q){p 2 — pq + q 2 ) = 0, 

(q 2 — p 2 )(uu + uu) + q(p 2 — l)(uu + uu) — p(q 2 — l){uu + uu) 
-5 2 (p 2 -q 2 )(p 2 -l)(q 2 -l)/(4pq)=0, 

(h(p)f(q) - h(q)f(p))[(uuuu + l)f(p)f(q) - (uu + uu)] 
+ (/b) 2 /(g) 2 - l)[(uu + u%)f(p) - (uu + m)f{q)} = 0, 



(Q2) 
(Q3) 
(Q4) 



where h 2 = f A + 5f 2 + 1. This can be parameterized with Jacobi elliptic functions: 
f(x) = k± sn(x, k), h(x) = sn'(x, k), 5 = —(k + 1/k). The simpler form of Q4 given above 
was actually discovered later in [5]. 

The A-series is auxiliary in the sense that Al goes to Ql with u — > (— l) n+m u and A2 
to Q3s=o by it — > w(~ 1 ) n+m . We will not discuss them here. 

2.2 Hirota's bilinear method 

Our approach to soliton solutions is based on Hirota's bilinear method[6j. It has been very 
successful in the continuous case and it is expected to be equally fruitful in constructing 
discrete multi-soliton solutions. Hirota's idea was to make a dependent variable transform 
into new variables, for which the soliton solution would be given by a polynomial of 
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exponentials. In terms of these new dependent variables the dynamical equations were 
quadratic and derivatives appeared only in terms of Hirota's bilinear derivatives. 

But what would be the natural discrete generalization? The essential property of an 
equation in Hirota's bilinear form seems to be its gauge invariance (cf., (?] in the continuous 
case) and it has a natural discrete extension, leading us to the following definition: 

Definition 1. We say an equation is in Hirota bilinear (HB) form if it can be written 
as 

^2 °i + u j"' m + f4) 9j( n + + Hj) = (2.2) 

3 

where the index sums z>y + vj = v s , ^ + [ij = fi s do not depend on j. (The functions 
f, g may be the same.) 

Proposition 1. Equations in HB form are gauge invariant, i.e., if functions fj,gj solve 
a set of equations in HB form, then so do the gauge transformed functions 

/j(n,m) = A n B m f 3 {n 1 m) 1 g'^m) = A n B m 9j (n,m). (2.3) 

Proof. We find 

f' 3 (n+v+ } m+^) g'^n+u^m+fxj) = A 2n+U t +V i B 2m+ ^ + ^ f^n+vf, m+/i+) gj {n+vJ , m+fij) 

but since the overall factor is the same in each term of the j-sum in (12. 2ft it can be taken 
out. □ 

For integrable equations in HB form there is a perturbative technique which leads to 
multi-soliton solutions, more or less algorithmically. This is described in the next section. 



2.3 Constructing background solutions 

First we have to construct the background or vacuum or seed solution, on top of which the 
soliton solutions are constructed. To do this we use the fixed-point idea [3] in which the 
consistent equations (12.11) are used with the assumption that u — u. However, since some 
equations are invariant under u — > T(u) it is actually sufficient that u = T(u). We only 
consider global invariances (that is, the transformation T is independent of n, m) and in 
order to keep the multi-linearity we assume that the transformation is linear fractional: 

u nm -> T(u nm ) := ° lUnm + ° 2 ; Cl c 4 - c 2 c 3 ^ 0. (2.4) 

C?,U nrn + C 4 

Depending on the equation there will be conditions on the parameters q. It is straight- 
forward to find the invariances of the equations, they are given in Table [U Note that the 
special cases with 5 = have a bigger invariance group. For Q4 there are various special 
cases depending on which limiting case of the elliptic parameterization one chooses. 
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Eq. 


Tin) 


T{u) when 6 = 


HI 


u + c, —u + c 


NA 


H2 


It 


NA 


H3 


u, —u 


cu, c/u 


Al 


u, —u 


cu, c/u 


A2 


u, —u, 1/u, —1/u 


NA 


Ql 


u + c, —u + c 


full Mobius 


Q2 


u 


NA 


Q3 


u, —u 


cu, c/u 


Q4 


u, —u, 1/u, —1/u 


various 



Table 1: Invariances ( 12.4}) of the equations in the ABS list. 



For each invariance T of a given equation we then need to solve 

Q(u,u,T(u),T(u);p,r) = 0, (2.5a) 
Q(u,u,T(u),T(u);q,r) = 0, (2.5b) 

in order to obtain the corresponding background solution (here r is a parameter of the 
solution). Such a solution then automatically solves (12. lap by virtue of CAC. 

It may be possible to construct still further solutions that can be called "background 
solutions" , but we will not consider them here. 



2.4 Constructing 1-soliton solutions 

Once the background solution is obtained we construct a one-soliton solution (1SS) using 
the CAC cube once more, with u now being the 1SS. This amounts, in fact, to a Backlund 
transformation (BT). 

Thus the first task is to solve 

Q(u,u,u,u;p,r) = 0, (2.6a) 

Q{u,u,u,u] q,r) = 0, (2.6b) 

where we take 

u = uo + v. (2.7) 

Here uq is the background solution and the bar-shift stands for possible modifications to 
get the f-equations into a suitable form. Then solving for v and v from (12.61) we get (for 
quadratic Q) 

- EV ° V (2.8) 



v + F 1 v + H 7 

where E, F, G, H may depend on n, m. The bar-shift modification in u$ is to be chosen 
so that the numerators have no constant term. By introducing v = g/ f and $ = (g, f) T 
we can write (12.81) as a matrix equation 

$(n + 1, m) = Af(n, m)§(n, m), $(n, m + 1) = M(n, m)$(n, m), (2.9) 
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where 

A/-(n,m)=A(f M(n,m)=A / ^ , (2.10) 

where E, F, G, H may depend on n, m. The constants of separation A, A' are to be chosen 
so that $ = $. In all cases studied in this paper it turns out that 

U a' m A, M(n,m)=i \, (2.11) 

where S, T, a, r, A, Q are constants and U nm some function of n, m (sometimes U nm = 1). 
If we now introduce ^ by 

*(n,m) = ( U * m f)*(n,m), (2.12) 



then for \l/ we have 



fs o\ (t o\ 

tf(ra + l,m)=l A J^(n,m), #(n,m+l)=( ^l*(n,m). (2.13) 



These are compatible if 



S-A T-Q' 

and in that case it is easy to derive 



From this we can construct m) and then v: Introduce 



and then v can be written as 



Now redefining the constant p ,o we can also write this as 

Pn,r 
1 + p 



(2.14) 



(S n T m \ 

T T n (s n T m A n n m ) A n n m ) *(°'°)- ( 2 - 15 ) 



= ' 4) (v) ) /; "- u (2 - 1G) 



u n,m U0:0 Pn,m/ P0,0 
^n,m — . r fo.o / / iV 



'n,m — U nm ■ (2.18) 
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2.5 N-soliton solutions and Casoratians 

Having a 1SS in the form mentioned above allows us to propose a change of dependent 
variables such that the original equation is given in terms of discrete equations in HB form. 
These equations will then be shown to have solutions given in Casorati determinant form, 
corresponding to the Wronskian form solutions for continuous HB equations. We will now 
discuss some generalities about the Casoratians. 

In general the Casorati matrix is constructed as follows: given functions ipi{n, m, I) we 
define the column vector 

^(n, m, I) = (ipi(n, m, /), 7/> 2 (n, m, I), ■ ■ ■ , ^(n, m, l)) T , (2.19a) 

and then the generic N x N Casorati matrix is composed of such columns with different 
shifts li, with the determinant 

Cn.mO; {k}) = |V>(n,m,Zi),V(fi,m,Z 2 ), • ■ ■ ,ip{n,m,l N )\. (2.19b) 
Two typical Casoratians that will be used later are 

c l,mW ■= \1>(n,m,0),il>(n,m, 1), ■ • ■ ,ip(n,m, N - 1)| 

= |0, 1, - - - ,N-1\ = \N~^1\, (2.19c) 

C 2 n ,mW := |V(n,m,0),--- ,il>(n,m,N-2),il>(n,m,N)\ 

= |0, 1, — , iV — 2, JV| = \N^2,N\, (2.19d) 

where we have also introduced the standard shorthand notation [9] in which only the shifts 
are given, and where furthermore sequential changes in the column index are indicated 
by a hat: |0, • • • , M, . . . | = \M, . . . |. We also denote \N^1\ = |1, 2, • • • , TV — 1 1 . 

In the Casoratians used in this work the entries ipi in the iVth-order vector (12. 19aj) are 
given by 

Vi(n, m, I) = g+(a + k t ) n (b + ^) m (c + k,) 1 + g r(a - h) n {b - ^) m (c - k t ) 1 , (2.20) 

where gf and ki are parameters (and in some cases c = 0). In this form the shifts in I 
(bar-shifts) are not in any way special and therefore we can equally well define Casoratians 
where the shifts are in n (tilde-shifts) or in m (hat-shifts). For later use we indicate these 
three different shifts by the shift operators E u , v = 1,2,3, respectively, i.e., 

E x ip = E 2 ^ = J, E^i) = ip. 

Down shifts are denoted by E u , v = 1, 2, 3, respectively. 

We can now define a Casoratian w.r.t. the three different kinds of shifts, 

| jV^l| M = E^, (^)V, • • • , (^)^VI, (y = 1, 2, 3). (2.21) 

Since 

(a M - a v )ip = (E^ 1 - E v )ip, /i,z/ = l,2,3, (2.22) 
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where 

«x = a, a 2 = 6, a 3 = c, (2.23) 

we have 

(E u ) j vp = [E» - (a„ - or„)]ty, fi,u=l, 2, 3, 
and substituting this into (12. 211) yields 

\N^l\ m = \N^\ m = \N^l\ m - (2-24) 

When the Casoratians are constructed using ip of (12.201) the size on the matrix iV 
indicates the number of solitons and the set {ki\f =l provides the "velocity" parameters of 
the solitons, while the parameters {pf,p^}i are related to the locations of the solitons 
(by gauge invariance only their ratio is significant). 

We shall here mention that the bilinear equations we get in the paper are more or less 
similar to the Hirota-Miwa equation [Til [T2] 



a(b — C)T nim j + iT n+ i t m + ij + b{C — a)T ni m+i t iTn+l,m,l+l + c { a ~ &)Tn+l,m,j7n,m+l,J+l — 0, (2.25) 

(cf. (15. 17115.201) ). or belong to the bilinear equations (13] which were derived by imposing 
the transformation [T2] 

oo j 
8=1 ■* 

on Sato's bilinear identity. The above transformation, referred to as Miwa transformation, 
provides a connection between continuous coordinates {xj} and discrete ones {lj}, and 
transforms the basic continuous plane wave factor 



exp 

into the discrete one 



oo 



x jPi + z2 x jQi (2.27) 



jx=l j=l 



OO 



CLjPi 



(2.28) 



The discrete exponential function (12.281) corresponding to the plane wave factor pi (1-2.2), 
plays a central role in the discrete r function in Hirota's exponential-polynomial form, 
while in Casoratians its counterpart has the form (I2.20p . 

In most cases a bilinear equation that can be solved by an N x iV Casoratian (or 
Casoratians) is reduced to a Laplace expansion of a 2N x 2N determinant with zero 
value. For later convenience we give the following lemmas. 

Lemma 1. |2J/ 

N , N v 

2j |ai, ■ ■ • ,a,-_i, ba i; a i+i , • ■ ■ ,a*| = (5> ) |ai, ■ • • ,& N \, (2.29) 

where a^ = (ay, ■ ■ • , aNj) T and b = (b%, ■ ■ ■ , &at) t are N -order column vectors, and ba^ 
stands for (6iay, • ■ ■ , bNaNj) T ■ 
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Lemma 2. [9] Suppose that B is an N x (N — 2) matrix and a, b, c, d are iV '-order 
column vectors, then 

|B,a,b||B,c,d| - |B, a, c||B, b, d| + |B, a, d||B, b, c| = 0. (2.30) 

In fact, the l.h.s. of (I2.30p is just the Laplace expansion of the following 2N x 2N 
determinant, 

1 

2 

3 HI 

3.1 Background solution and 1-soliton solution 
3.1.1 The background solution 

The HI equation is given by 

HI = (u — u){u — u) — (p — q) = 0. (3.1) 

Using the fixed point idea with transformation u = T(u) = u + c we get the side-equations 
of the CAC-cube in the form 

(u — u) 2 = r + c 2 — p, (u — u) 2 = r + c 2 — q. 

For convenience we absorb r into c 2 , and reparameterize (p, q) — > (a, b) by 

p = c 2 - a 2 , q = c 2 - b 2 . (3.2) 

and then the above equations factorize as 

(u — u — a) (u — u + a) = 0, (u — u — b) (u — u + b) = 0, (3.3) 

Since the factor that vanishes may depend on n, m we actually have to solve 

u-u=(-l) a, u-u = (-l) x b } (3.4) 

where 6 1 , x G Z may depend on n, m. Furthermore, since the value of the exponent is 
only relevant modulo 2 and n 2 = n mod 2, the exponents must be linear combinations 
of 1, n, m, nm with coefficients or 1. 

The integrability condition for f)3.4p leads to 

b[(-iy-(-iy]=a[(-i) e -(-if], 

and since a, b are independent this just means that o = cr(m), 6 = 9(n). Since these 
are defined modulo 2 we can have 9 = tin + t , a = Sim + s where tj,Sj G {0,1}. 
Furthermore, since a, b were only defined up to sign we may take t = s = 0. 
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B a b c d 
B a b c d 



There are now two essentially different cases for each exponent, either 6{n) = or n, 
and a(m) = or m. The choice 8 = leads to a solution u = an + . . . while 9 = n leads 
to u = — 1(— l) n a + . . . . These solutions can be combined, and thus we get the following 
set of possibilities for u° 

an + bm + 7, (3.5a) 
|(-l) n a + 6w + 7, (3-5b) 
an+|(-l) m 6 + 7 , (3.5c) 
\{-l) n a + \{-l) m b + 1 . (3.5d) 

We could also try to get a background solution using the fixed point idea with the 
other transformation mentioned for HI in Table [H namely u = T(u) = —u + c, but it 
produces the same set: with a reparameterization 

p = r + a 2 , q = r + /3 2 , u = y + |c. (3.6) 

the side-equations become 

(y-y- a){y-y + a) = 0, (y-y- P)(y-y + (3) = 0, (3.7) 

which is as in (13. 3p . The constant c can be absorbed into the constant 7 and the free 
parameter r redefined when comparing fl3.2IUI3.6l) . and thus the background solutions are 
as in (13.51) . 

It seems that only (I3.5ap leads to a soliton type solution so in the following we only 
consider it. 

3.1.2 Constructing 1SS using a BT 

The BT generating 1SS for HI is 

(u — u)(u — u) = p — x, (3.8a) 

{u — u)(u — u) = x — q. (3.8b) 

Here u is the seed solution f)3.5ap . x is the parameter in the bar-direction, and we search 
for a new solution u of the form 

u = u + v , (3.9) 
where uq is the bar-shifted background solution (I3.5al) : 

u = an + bm + k + 7, (3.10) 

where k is related to x by 

X = C 2 _A; 2 . (3.11) 
Substituting fl33al) . fl3H and flBTTUj) into (jMj) yields the desired form 

Ev ^ Gv . n , 

v + F v + H 
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where 

E = -(a + k), F = -{a-k), G = ~(b + k), H = -(b-k). 

Now matrices A/", M. defined in ( I2.10p are n, m independent and since E — F = G—H = 
—2k they commute and we can take A = A' = 1. Then, following the method presented 
in Section l2~4l we find 

g!g^rg F n H m J $(0,0), (3.13) 



and if we let 

'E\ n (G\ m fa + k\ n fb + k 



then we obtain 



Pn,m=(j Po,o= — H ^) Po,o, (3.14) 



fit 



-2kp r 



1 + Pn,' 



(3.15) 



Finally we obtain the 1SS for HI: 



«S = (an + 6m + 7 ) + A; + - f P "' m (3.16) 



3.2 Multi-soliton solutions 

In an explicit form the 1SS given above [ (I3.16P with (I3.14p ] is 



1 H~~ pn , m 

fc[e+(q+fc)"(b+fc) m -g-(a-fc)"(b-fc) r 



an + 6m + 7 — 



(3.17) 



e+(a+fc) n (M-fc) m +e-(a-£i) n (6-fc) m ' 



where we have separated po,o to q~/q + . 

The above suggests that the basic ingredient in constructing an NSS is ip of (12.201) 
with c = 0, i.e., 

^{n, m, I- h) = gfkl(a + h) n {b + k t ) m + gr(~k t ) l (a - h) n {b - h) m . (3.18) 

Here the index I is needed as a column index in the Casorati matrix. 

The form of the 1SS (13 . 1 7[) suggests a generalization to NSS, and after checking ex- 
plicitly the 2SS and 3SS, we arrive at the following proposition: 

Proposition 2. N-soliton solution for HI is given by 

<» = «. + to, + 7 -£, (3.19) 

where f — \N — 1| [3] , g = \N — 2, N\ , with ip given by ( I3.18p . 
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Note that this solution is the same as (5.29) of Part I with A = 0. Indeed using 
(1-2. 1,2.3,2.8c) one finds that 

s (0 - 0) -l> = i 

and the parameters gf in (13. 18j) and Cj in (1-2.3) are related by 

^ = ^_J]^L^. (3.20) 

£?,- 2fcj kj -\- ki 

However, the solution (13.191) contains more freedom in the form of the additional param- 
eter c, cf. (|3I2D and (1-5. Id). 

3.2.1 The bilinear form of HI 

In order to prove Proposition [2] we write (13.11) in bilinear form. First we introduce a 
dependent variable transformation 

<l 8 = an + bm + j- 9 f^. (3.21) 

Jn,m 

When this is substituted into (13. ip we get a rational expression quartic in /, g. In order to 
split this expression into bilinear equations we note that if / = \N — 1| , g — \N — 2, N\. , 
then for N = 2 one can quickly scan for possible equations solved by them from the set 

aogf + a{gf + a 2 gf + a 3 gf + a A fJ+ a 5 ff + a 6 gg + a 7 gg = 0, 
and find the bilinear equations 

H x = gf-gf + {a-b)(ff-ff) = 0, (3.22a) 
H 2 = <7/-f/+(a + 6)(//-/7) = 0. (3.22b) 
After this it is easy to show that 

HI = ( u - - u) - p + q = - [Hx + (a - b)fj] [H 2 + (a + b)ff] /(ffff) + (a 2 - b 2 ), 
and thus the pair (13.221) can be taken as the bilinear form of (13. ip . 

3.2.2 Proof of the Proposition d 

Now it remains to show that the /, g given in Proposition [2] solve equations (I3.22p . We 
prove (I3.22ap in its tilde-hat-downshifted version, 

9f-9f+(a-b)(ff-ff) = 0. (3.23) 
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Comparing ipi(n,m,l) given by (13.181) and (lA.lj) . we can use the formulas given in Ap- 
pendix |A] with c = 0, and in this subsection we always have shifts in the /-index, i.e., we 
take k = 3 for the formulas in Appendix [A] Let us write (13.231) as 

- (a - b)ff + f(g + af) - f(g + bf) = 0. (3.24) 



In this formula, / = \N — 1| [3] , and for f,f,g + af, f and g + bf we use ( lA.6kl) with 



\i = 1 and v = 2, (1A.6al) with /i = 2, (1A.7al) with /i = 1, (1A.6al) with /i = 1 and (lA.7al) 
with /i = 2, respectively. Then we have 

a N-2 b N-2 { _ {a _ b)f f + f { g + _ ^ + hf )] 

= -\^l\ [3] \N^3,iP(N -2)^(N -2)\ [3] 

+ \N^2, f(N - 2)| [3] \N^3, N - 1, f(N - 2)| [3] 
-|A^2,^(iV - 2)| [3] |iV^3,iV - l,^(iV - 2) l[3j 



0. 



where we have made use of Lemma [2] in which B = (N — 3), and a = ip(N — 2), b = 
■0(iV - 1), c = f(N - 2) and d = f(N - 2). 

For ( I3.22bl) we use its tilde-down-shifted version 

9j- 9f + (a + &)(//- ff) = 0, (3.25) 

and rewrite it as 

-(a + b)ff + f(g + af) - f(g-bf) = 0, 

where / = \N — 1L, and for f,f,g + af, f and g — bf we use ( IA.61I) with /i = 2 and 
v = 1, (IA.6cD with /i = 2, HA.7a|) with /i = 1~ flA.6al) with /j, = 1 and (IA.7bp with /i = 2, 
respectively. Then we have 

' N-2rN-2 



ar-W-\-lp + b)ff + /( 5 + af) - f(g-bf)] 



|0 2 | 

-|jV-l| M |jV-3,^(JV - 2), EV(iV - 2)| [3] 
-|jV-2, E^{N - 2)| [3] |jV^3, iV - 1, f(N - 2)| [3] 
+ | iV~^2, ^>(JV - 2)| [3] | N^3, N - 1, E 2 ^(N - 2 
0, 



where we have made use of Lemma [2] in which B = (N — 3), and a = ip(N — 2), b = 
ip(N-l), c = i)(N-2) a,ndd = E 2 i)(N -2). □ 
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4 H2 

4.1 Background solution 

The equation H2 is given by 

H2 = (u — u)(u — u) — (p — q)(u + u + u + u + p + q) = 0. (4.1) 
After reparameterization 

p = r — a 2 , q = r — b 2 , (4.2) 

the equations on sides become 

(u - uf - 2a 2 (u + u)+ a 2 (a 2 - 2r) =0, (u - u) 2 - 2b 2 (u + u) + b 2 (b 2 - 2r) = 0. (4.3) 
After the further substitution 



u = y 2 - \r (4.4) 



the equations ( 14.3}) factorize as 



(y + y + a){y + y - a)(y - y + a)(y - y - a) = 0, (4.5a) 
(y + y + b)(y + y-b)(y-y + b)(y-y-b) = 0. (4.5b) 

Thus we have to solve 

y-(-l)*y = (-l)» 0> y-(-l)'y = (-l)*6. (4.6) 

where the exponents a, 6, p, <ft G Z are again some linear combination of 1, n, m, nm with 
coefficients or 1. Equations (14.61) are satisfied if 

p = ns 2 +r 1 , a = ms 2 + si, 9 = m p+a+n(t n +si)+ti, (j) = na+p+m(p m +r 1 )+p 1 , 

but most of the freedom is superfluous: Since the sign of y was left undetermined in (14.41) 
we can redefine 

„. . ( _-\\nm S2+n si+mri 

yn,m v 1 ) yn,m 

and then the equations simplify to 

y-y = (-ir'" +tl a, y-y= (-1)"»^^6, 

already analyzed after (13. 4p . Thus we find the possible background solutions u^fm of the 
form 

[an + bm + 7] 2 — |r, (4.7a) 
[\{-l) n a + mb + 1 ] 2 -±r, (4.7b) 
[na + |(-l) m 6 + 7 ] 2 -ir, (4.7c) 
[i(-ira+i(-ir6 + 7 ] 2 -|r. (4.7d) 
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4.2 1-soliton solution 

Note that H2 is invariant under simultaneous translation of p, q by t and u by — \t. We 
use this freedom to eliminate r in order to simplify the presentation. 

Now we take (j4.7a|) as the seed solution to construct for H2 through its Backlund 
transformation 

(u — u)(u — u) = (p — x)(u + u + u + u + p + x) , (4.8a) 
(u — u) (u — u) = (x — q)(u + u + u + u + x + q) . (4.8b) 
where u is the seed solution (I4.7al) and we search for a new solution u of the form 

u = u Q + v 1 (4.9) 

where uq is the bar-shifted background solution (I4.7al) : 

u = (an + bm + k + >y) 2 , (4.10) 

with 

K=k 2 . (4.11) 
Substituting these into (gSJ yields fl2TT0|) with 

E = -2(k + a)[(n + I) a + 777,6 + 7], F = -2(k - a)[na + mb + 7], 
G = -2(k + 6)[na + (m + 1)6 + 7], H — -2{k - b)[na + mb + 7]. 

and the corresponding matrices J\f, Ai are compatible if we take 

A = A' = —l/(2U nm ), U n>m := an + 6m + A. (4.12) 

With this U we obtain (12X0 with 

5 = a + A;, A = a - k, T = b + k, Q = b-k, a = t = -1/2. (4.13) 

Then defining 

fS\ n fT\ m f a + kYfb + k\ m 

p ^={a) UJ Po ' = l^J IrntJ Po ' ' (414) 

we find 

4:kU ntm p n ^ m 



1 + /V, 

and finally we get the 1-soliton for H2 in the form 



(4.15) 



<L 5 = + 2^ n , m i-^ + k\ (4.16) 

-L ~r Pn,m 



15 



4.3 Multi-soliton solution 

Motivated by the structure of 1SS (14.161) (cf. (13.171) ). and after checking 2- and 3-soliton 
solutions we propose the following Casoratian expression for the NSS: 

2 \N^~2,N\ m \N^3,N-l,N\ m + \N^~2,N + l\ m 

\N — 1 AT — 1 

l JV ^[3] l JV X l[3] 

where the matrix entries are as for HI ( 13.181) . In order to prove this we derive a bilinear 
form of H2. We propose 



NSS_TfZ _r > jj (A-\R\ 
u n,m ~ u n,m iu n,m j. 1" - , 

where /, h and s should satisfy 

h-s = >yf, (4.19) 

where 7 is some constant. Indeed, if 
f=\N^l\ w g=\N^2,N\ [3] , 3=\N^3,N-1,N\ W h=\N^2,N+l\ w (4.20) 

then noting that (E 3 ) 2 ipi(n, m, I) = ip^n, m, I) = ipdn, m, I + 2) = kft^^n, m, I) and using 
Lemma [H we have (J2iLi kf)f = h — s, i.e. 7 = J2i=i tf- 

The solution (I4.18P is the same as (1-5.26) with A = and 

s m = J2jkj + g - f 2SW - 2(£ i ^ (0 ' 0) + (£^) 2 = ^p, 

which can be found by using (1-2. 1,2. 2, 2. 3, 2. 8c, 5. 23). 

Under condition (I4.19p . H2 can be represented through the following bilinear system, 

Hi = gf-gf + (a-b)(ff-fj) = 0, (4.21a) 

H 2 = g?-gf+(a + b)(ff-ff) = 0, (4.21b) 

H 3 = -(a + b)fg + a% + bft + }h-fh = 0, (4.21c) 

U 4 = _( a _ b ) f ~ +a f§_ b fg + f h-fh = 0, (4.21d) 

H 5 = b{fg-fg) + fh + fs-gg = 0, (4.21e) 

in which 7i\ and 7^2 already appeared in (13.221) . In terms of the above bilinear equations 
H2 can be given as 



H2 = Y^HiPi, (4.22) 
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with 



Pi 

P2 
P3 
P4 
P, 



-4(a + b) (UU - a 2 + b 2 )ff - Ufg - (a + b)fg 



a _ h )(pu _ a 2 + 6 2 } JJ + ( ££ _ fl2 + &2) J ? _ _ (a _ b p f ~ 

(a - b)Uff + Ufg - -fh + fh 
(a + b)(Uf7-fg) + U(7g-m 



4(a 2 - 6 2 )//, 

where U was defined in (14.121) . It remains to prove the following 

Proposition 3. The Casoratian type determinants f , g, h and s given in (14.201) with 
entries given by ( 13.181) solve the set of bilinear equations ( 14.211) . 

Proof. Among (I4.2ip TCi and Ti.2 have been proven before. Next we prove (I4.21cl) in the 
following form 

- (a + b)fg + f(h + ag) - f(h - bg) = 0, (4.23) 

which is a down-tilde-shifted version of the original one. Since ipi(n,m,l) given by ( 13.181) 
is just (IA.1I) with c = 0, we use the formulas given in Appendix [A] with c = and k = 3. 
In ( Km g = \N^2, N\ [3] , and for f,f,h + ag, f and h - bg we use (TOIh with /i = 2 
and v = 1, (IA.6cp with \i = 2, (IA.8aP with /i = 1, (IA.6aD with /i = 1 and (IA.8bD with 
fi = 2, respectively. Then we have 

1 „N-2rN-2 



1^2 



— |iV — 2, AT| [3] I jV - 3, V(iV - 2), £ 2 ^(iV - 2) 



[3] 



-\N - 2, £ 2 </>(N - 2)| |JV - 3, N, tp(N - 2 



I [3] 



+\N - 2, f(N - 2)| [3] |JV - 3, N, E 2 ^{N - 2)| 
0, 



[3] 



where we have made use of Lemma [2] in which B = (N — 3), a = ip(N — 2), b = ip(N), 
c = f(N - 2) and d = E 2 ^(N - 2) 

( j4.21dl) can be proved similarly after a down-tilde-hat-shift. 

Next we prove (14. 2 lei) , which can be written as 



f(h-bg)-g(g-bf) + fs = 0, 



(4.24) 



where / = |JV- 1| [3]) g = \N -2,N\ m , s = \N — 3, N — l,N\ m , and h - bg, g-bf and 



I [3] 



[3] ' 



/ will be provided by (IA.8bl) with fj, = 2, (IA.7bj) with /j, = 2 and (IA.6cP with fj, = 2, 
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respectively. Then we have 

1 



b^y(h-bg)-g(g-bf) + fs] 



\n 2 \ 

\N^l\ [3] \N^S 1 N,E 2 ^(N-2)\ [3] 
-\N^~2,N\ [3] \N^~3,N-l,E^(N-2)\ [3] 
+ \N^2, E 2 ^(N - 2) | [3] |iV^3, N — 1, N\ [3] 



0. 



where use has been made of Lemma [2] with B = (JV — 3), a = ip(N — 2), b = ip[N — 1), 
c = ^(N) and d = £ 2 ^(iV-2). □ 



5 H3 

5.1 Background solution 

H3 S is given by 

H3 5 = p{uu + uu) — q{uu + uu) — 5(q 2 — p 2 ) = 0. (5-1) 
The side equations for T(x) = x then read 

r{u 2 + u 2 ) — 2puu = 5{p 2 — r 2 ), r{u 2 + u 2 ) — 2quu = 5{q 2 — r 2 ). (5-2) 
In this case we reparameterize 
p = r cosh(a') , q = r cosh(/3') , u nm = Ae Vnm + Be~ Vnm , AB = -\r5 (5.3) 
and then the equations (15. 2p factorize as 

( e y- y+a ' - l)(e^ y - a ' - l)( e y +y+a '- hl % - l)(e y+y - a '- ln % - 1) = 0, (5.4a) 
( e ^-y+^_ 1 )( e y-y-^'_ 1 )( e ^+/3'-in|_ 1 )( e ^-/3'-i n |_ 1 ) = Q ^ 5Ah ^ 

Since we only consider real u the various possibilities can be represented as in (14.61) . The 
analysis is then the same, especially since also here the sign of y is undetermined in (15.31) . 
Thus the solution for y is as in (13. 5p and for u we have 

u 0SS _ A e Vnm _|_ J^ e -y nm _ ^ e a'n+/3'm+7 _|_ j^^-a 1 n-fl 

= Aa n (3 m + Ba~ n (3- m , AB = -\r6, (5.5) 

where a = e a ' , j3 = e@' . 

In the case of T(x) = —x, one can find that the side equations are just as for T(x) = x 
except that r — > — r. Since r is a free parameter this adds nothing new. 
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5.2 1-soliton solution for H3 S 

Now we take (I5.5P as the seed solution to construct for H3 5 through its Backlund 
transformation 

p(uu + uu) — h{uu + mm) = 8(x 2 — p 2 ), (5.6a) 
x(mm + uu) — q(uu + uu) = 5{q 2 — x 2 ), (5.6b) 
where u is the seed solution (15. 5ft . and we search for the 1SS u of the form 

u = u + v 1 (5.7) 

where uq is the bar-shifted background solution (15.51) : 

u e = Aa n (3 m K + Ba- n (3- m K-\ (5.8) 

with x and k related by 

1 + K 2 . . 

k = t^—. (5.9) 

Following again the procedure in Section 12.41 we get JV, Ai in the form (12. lip with 

1 — a 2 /t 2 —a 2 + k 2 
S = r— , A = r — , a = p, 

ZOLK ZCtK 

1 - (3 2 k 2 -/? 2 + /t 2 



T = r , Q = r , r = q, 

U n , m = Aa n (3 m - BoT n fi- m . 



Then defining p by 



S y/ T \m / a V-l\7A 2 -, 
*"»=(aJ UJ ^U^V Po ' ' (5 ' 10) 



we find 



and finally 



,mPr, 



1SS Aa"/3 m (l + K ~ 2 p n , m ) + BoT^/T^l + K 2 Pn , m 



(5.11) 



J- ~r pn,m 

5.3 Bilinear form and Casoratian solutions 
5.3.1 iV-soliton solution 

Noting that p n ^ m given by ( I5.10P is in a "twisted" form in comparison with ( I3.14p and 
(14. 141) . we first introduce the Mobius transformations for the parameters 

2 a ~ C p2 b — C 2 k — C 2 r2 ° 2 2 r2 ° 2 /r 10 \ 

a =-—c^ = -bT~c K =-kT-c P= c^> q= c^> (5 ' 12) 
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which also contain a new auxiliary parameter c. This brings p n ^ m of (I5.10P into the 
canonical form 

' a + k\ n ( b + k s 



Pn, m =[^) [^) Po,o- (5.13) 
In terms of the new form of p n ^ m we can write the 1SS (15. lip as 

where ip is as in (12.201) . 

On the basis of the above 1SS, we propose that the NSS of H3" 5 can be given by 

7 f 

'• vs ' s - Aa n (3 mJ - + Ba~ n p- m y, AB = -\r5, (5.15) 



where / = \N — lL, with entries (12.201) . We may consider (15.151) as a dependent variable 
transformation for (15.11) . 

The solution (15.151) is the same as (1-5.21) with B = C = 0. In fact, let - in ( 15.151) 
equal to a which is the direction parameter for the bar-shift in Part-I. Then comparing 
(1532]) with (I-5.1c), using (1-2.32,5.19) and substituting p i in (1-2.2) by 

( p + k i \ n f q + h \ m ( o + fcA ' ^ lg . 

\P-ki) \q-h) \a-ki) l ' 



one finds 



^ = «- n /3" m , no)/n i (a-*i) = 7, V(-o)x 



j, V(-o)xn i (a-*i) = j: 



with the same parameter identification as in (13.201) . 

5.3.2 Bilinearization-I 

After introducing the two bilinear equations 

Bt = 2cff + (a-c)Jf-(a + c)Jf = 0, (5.17a) 
B 2 = 2cff+(b-c)Jf-(b + c)Jf = 0, (5.17b) 

we can represent H3 5 ( 15.11) as 



H3 



, _ -a 4 " +2 /3 4m+2 (a + c)(6 + c)<5 2 Pi + 4a 2n /? 2m <5P 2 P 2 + 16(a + c)(6 + c)P 4 P 3 



32a 2 «+ 2 /3 2m + 2 (a + c) 2 (6 + c) 2 B 2 ffff 
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where 



Pi 
P2 
P3 



f 
2c 

/ 



(b-c)fB 1 -{a-c)fB 2 



f 



(b + c)fB 1 -(a + c)fB 2 



+ c)(6 - c)(ffB 1 + ffBx) - {a + c)(a - c)(ffB 2 + ffB 2 ) 



(b + c)fB 1 - (a + c)fB 2 



f 



(b- c )fB 1 -{a-c)fB 2 



Thus (15.171) can be considered as a bilinearization of H3 5 , and the final step in constructing 
the NSS is 



Proposition 4. The Casoratian 

/=|JV^1| H , (v = 1,2 or 3), 

wii/i entries given by ip of ([27201) . sofoes tfie WKraear i/5 5 (157171) 
Proof. We prove (I5.17al) in its down-tilde-shifted version 

2cff + (a-c)Jf-(a + c)Jf = 0. 



(5.1* 



(5.19) 



For this equation we use Casoratians w.r.t. hat shift, i.e., k = 2 in (1A.6|) . /, /, /, /, / 
and / are given by the formulas (!A.6aP with // = 1, (1A.61j) with fi = v = 3, (!A.6cP with 
/i = 3, (1A.6ki) with /i = 1 and v = 3, (1A.61|) with /i = 3 and z/ = 1, and (!A.6aP with /i = 3, 
respectively. Then we have 

1 ( a _ 6) "-2 (c + 6) ^-2 (c _ b ) N - 2 [2cff + (a - c)// - (a + c)Jf] 



|r 3 | 

-|iV^2, V(iV - 2)L |JV - 3, ^(iV - 2), #ty(W - 2)| 



[->] 



+|JV - 2, E 3 ^(N - 2)| [2] |JV - 3, m - 2), ^(iV - 2)| [2] 

+ |A^2, </>(iV - 2)| [2] |iV^3, ^>(iV - 2), ^(iV - 2)| 
0, 



where we have made use of Lemma [2] in which B = (N — 3), and a = ip(N — 2), b = 
^p(N - 2), c = ip(N - 2) and d = E 3 ip(N - 2). 

The other bilinear equation ( 15. 1 7bh can be proved in its down-hat-shifted version in a 
similar way by taking k = 1. □ 

5.3.3 Bilinearization-II 

In fact there is another bilinearization of H3 5 using (15.151) . Consider the bilinear system 

B[ = (b + c)ff+(a-c)fj-(a + b)fj=0, (5.20a) 
B' 2 = ( c -b)7f-(a + c)fj + (a + b)fT=0, (5.20b) 
B' 3 = (c-a)(b + c)Jf + (a + c)(b-c)Jf + 2c(a-b)f7=0. (5.20c) 
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This system is related to H3 5 through 



m d 



ffff 



2 ^,2n g2m f fB[- f fB' 2 j_ _g2 Q -2n ^-2m /^~/ f&2 



A 01 {a+c)(b+c) 

ffB' 2 +fjB 2 TfBi+ffB^ 



+ ab( 



(a—c)(b—c) 
2(a+b)ffB' a 



(a+c)(b-c) (a-c)(b+c) (a 2 -c 2 )(6 2 -c 2 ) 



The bilinear system (15.201) shares the same Casoratian solutions H5. 1T[) : 
Proposition 5. The Casoratian 

f = \N^l\ W] , {u = 1,2 oi 3), 

with entries given by ip of fl2T2DD . solves the bilinear H3 5 (J52DD- 



(5.21) 



Proof. We only prove fl5T20aD and CTB . flpQbD is similar to flOUij) . We prove fl5T20aD 
in its down-tilde-shifted version, i.e., 



(b + c)ff + {a-c)ff-{a + b)ff = 0. 



(5.22) 



We need to use Casoratians w.r.t. bar shift. So we now fix k = 3 in (1A.6I) . For /, /, /, 
/ and /, we use the formulas flA.6dj) with /j = 2, ( A.6g ) with /j = 1, (IA.6bl) with fj, = 1, 

(IA.6ip with /i = 2 and (IA.6nO with // = 2 and u = 1, respectively, and / = |iV — 1| [3] . 
Then we have 



(a - c)"- a (6 + c)"- a [(6 + c)ff + (a - c )// - (a + b)ff] 



\N-2 



r 2 



■|iV-2,SV(JV-l)| [s] l^-l,^-l)l 



I [3] 



+|AT - 2, ^(iV - 1)| [3] | N - 1, E 2 tfj(N - 1 



I [3] 



-|JV - 1| |JV - 1, i>(N - 1), E 2 i>(N - 1)| 



[3] 



0. 



where we have made use of Lemma[2]in which B = (N — 2), and a = ip(0), b = if)(N— 1) 
c = f(N - 1) and d = E 2 ip(N - 1). 

For (I5.20cl) . after a down-hat shift we get 



2c(a - b)ff - (a - c)(b + c)ff + {b-c)(a + c)ff = 0. 



(5.23) 



In this case we fix k = 1 in (1A.6|) . For /, /, /, /, / and /, we use the formulas (lA.6bl) 
with /i = 2, (lA.6np with // = v = 3, (lA.6np with // = 3 and z/ = 2, (1A.6bl) with /i = 3, 
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(1A.6dj) with /i = 3, and (lA.6ml) with /i = 2 and v = 3, respectively. Then 
1 



(6 _ a f'\c + a)»-\c - a) N ~ 2 [2c{a - b)ff - (a - c)(b + c)ff + (b-c)(a + c)ff] 



|r 3 i 

= -|aTT2, f(N - l)\ m \N-2, 4>(N - 1), E 3 ?P(N - l)| [x] 
+ \N^2, i;(N - 1)| (1] |iV^2, ^(iV - 1), #ty(iV - 
+ |iV^2, E^(N - 1)\ {1] \N^2^{N - 1), f(N - l)| m 

= 0, 

where in Lemma [2] we take this time B = (N — 2), and a = ip(0), b = ip(N — 1), 
c = f(N-l) and d = £ 3 V>(iV- 1). □ 

6 Ql with linear background 

6.1 Background solution with T(x) = x + c 
Ql 5 is 

Ql 5 = p(w — u)(u — w) — g(w — «)(« — «) — 5 2 pq(q — p) = 0. (6.1) 
With the fixed point defined by T(x) = x + c the side equations for the background are 

r{u — w) 2 = p(c 2 + 5 2 r(p — r)), r(w — w) 2 = q(c 2 + 5 2 r(q — r)). (6.2) 

After the reparameterization (p, q) — > (a, b) with 

c 2 /r - 5 2 r c 2 /r - 5 2 r 

P= q2 _ s2 , g = b 2_$2 1 a:=pa, (3 : = qb, (6.3) 

equations (16.21) factorize as in (13.31) . and thus the OSS will be as in (I3.5p . where, however, 
we should replace a and b with a and /3, respectively. 

6.2 1-soliton solution 

The BT for constructing the 1SS is 

p(u — u) (u — u) — x[u — u) (u — u) = 5 2 p>c{x — p) ) (6.4a) 
k{u — u)(u — u) — q{u — u)(u — u) = 5 2 xq(q — x). (6.4b) 

Following the usual procedure we take u = an + (3m + 7 as the OSS and 

u = an + [3m + 7 + n + v , (6.5) 

as the 1SS, where v is to be determined. If we now choose 

c 2 /r — 5 2 r 



k 2 -5 2 
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k = kx, (6.6) 



then we get (12.101) with 

E = -x(a + k), F = —x(a — k), G = -x(b + k), H = — x(6 — k). 
Thus if we define 

a + k\ n fb + k 



Pnm -\a~k) \b^k ] '" ,n 

we find the 1SS in the form 

1 - p n 



u = an + (3m + 7 + k 



1 + Pr. 



where p, q, x depend on a, b, k as given in (16.3p . (16.6H . This is similar to ( 13.161) except for 
the more complicated dependence on the parameters a, b, k. 

6.3 NSS 

After studying the 2SS using Hirota's perturbative method we propose that the NSS is 
obtained using the following: 

= » + + ? -tf/r-^f (6-7) 

where 

/=|jV^i| M> 0=|-l,tf-l| M) (6.8) 

with ip defined by 

^(n, m, I) = g+(a + h) n (b + h) m {5 + h) 1 + g ;(a - h) n {b - h) m {5 - h)\ (6.9) 

The solution (16.71) is consistent with (1-5.11): Taking c = in (16.71) and A = D = 
0, B = 5/2 in (1-5.11), then using (I-2.34c,5.1b,5.4) one finds that 

a = 5, S(-a,a)='£ jl ± i + j. 

However, the additional parameter c (16.71) implies more freedom, cf. (16.31) and (1-5. lb). 

The functions /, g satisfy the following bilinear equations (see proposition below) 
among others 

Qi = ?f(b-S)+ff(a + S)-ff(a + b) = 0, (6.10a) 

Q 2 = Jf(a-b) + ff(b + 5)-ff(a + 5) = 0, (6.10b) 
Q 3 = -ff + fg(-a + 5) + ff + fg(b-6) + f${a-b)=0, (6.10c) 

Q 4 = %(a-b) + fg(a + b)-fg(a + b) + m-a + b)=0. (6.10d) 
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We note that Q4 can also be replaced by 

Q4 = Qs + Qi = If- Jg{a + 5) - J/ + 7£(5 + 6) + 7^(a - b) = 0, (6.11) 

which is similar to Q3. When the dependent variable transformation (16.71) is substituted 
into Ql 5 we find that the result can be expressed in terms of the Qi defined above: 

qi'= (c2/r ^ 2r)3 , fa,P„ (6.12) 

(a? - - P)(a - 6)(o + S)ff/ff i=i 

with 

Pi = (a-b)[ffg{-a + b)+f{fg-fg){a + b) 

+ fgg(-a 2 + 5 2 ) + jggib 2 - S 2 ) + fgg{a 2 - b% (6.13a) 
P2 = (a + b)[ffg(a-J) + f(fg~mb-S)+g(fg~fg)(a-b)(b-S)U6.l3b) 
p 3 = ( a + b )( a + S)tffg(a-b)j-ffg{b-6) + ffg(-a + fi)], (6.13c) 
Pi = f(a + 8)[ff(-a + b) + (fg-fg)(a-5)(b-8)\. (6.13d) 
Then it remains to prove the following: 

Proposition 6. The Casoratian type determinants f , g, given in (16.81) with entries given 
by (16.91) solve the set of bilinear equations (16 . lOj) . 

Lemma 3. By means of Lemma\^ the following formulas are zero, (fi = 1,2), 

Y,=f\E^(-l),-l,N^2\ [3] +£\E^{-l),I^l\ [3] -g\E^{-l),ff^2\ [3] , (6.14a) 
^3 =|^(-1), -1, iV^2| [3] |^(-1), iV^Ti[ [3] — |^(-1), -1, iV^=^| [3] i^^l| [3] 

+ g\f(-l),i ! (-l),N^2\ [3]] (6.14b) 

^ M =/|^V(— 1), -l,iV^2| [3] + /|^^^(-1),7V^1| [3] - t7|^V(-l),^V^2| [3] , (6.15a) 
Z 3 =|EV(-l),-l,^2| [3 J^(-l),iV^l| [3] 
- |EV(-1), -1, N^2\ [3] I^V(-l), N^l\ [3] 

+ g\i?M-l),E 1 4(-l),N^2\ [Si . (6.15b) 

Proof. For the above formulas, we can use Lemma [2] by respectively taking 
B = (iV^2), (a,b,c,d) = (^(-1),^(-1),^(0),^(JV- 1)), 
B = (N^2), (a,b,c,d) = (f(-l),f(-l)M-l)MN-l))-, 
B = (N^2), (a, b, c, d) = ^(-1), ^(0), ^{N - 1)), 

B = (N^2), (a,b,c,d) = (EV(-l),^V(-l),^(-l),^(iV- 1)). 

□ 
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Proof. 

Proof for (I6.10al) : The down-bar shifted (16.1 Oap is nothing but B' 2 of H3 5 with c = 5. 
Proof for ( I6.10bl) : By a down-tilde shift (I6.10bl) is written as 

(a - b)fj +(5 + 6)7/ - (a + 5)Jf = 0. (6.16) 

We fix k = 2 and c = 5 in (IA.60 . For /, /, /, / and /, we use the formulas (1A.6bl) with 
H — 1, ( 1A.6il) with /x = 3, ( lA.6dl) with /x = 3, ( A.6g ) with /x = 3 and (1A.6nl) with /x = 3 
and v — 1, respectively, and / = |iV — 1| [2] . Then we have 

' - 6)^ 2 (5 + 6)^- 2 [(a - 6)/7 + (* + 6)7? - (a + 5)7/] 



|r 3 r 

-|A^,i 3 ^-i)i [2] |iv^i^(iv-i)| [2] 

-\rf^i\ m \r^iMN-i),E 3 ip(N-i)\ l2] 
= o, 

where we have made use of Lemma [2] with 

B = (xV^2), (a, b, c, d) = (</>(0), ip(N - 1), ^(xV - 1), E^(N - 1)). 
Proof for (16. 10c[) : By a down-tilde-hat shift (16. 10c[) is written as 

- J[f + (a- 5)g] + /[/ + (6 - 5)g] + (a - 6)7<7 = 0. (6.17) 

This time we fix k = 3 and c = 5 in (1C.1I) . For /, /, g, f, f, g and /, we use the formulas 
flClbl with fi = 2, flClaj) with xx = 1, flOTcl) with /x = 1, flCTTB with /x = 1, flOTaj) with 
/x = 2, (IC.lcl) with /x = 2, and ( |C.lg ), respectively. Then we have 

-/[/ + (a - %] + /[/ + (6 - S)g] + (a - 6)/ 5 
= (a - 5) 2 F! - (6 - 5)V 2 + (a - <5) 2 (6 - 5) 2 Y 3 
with defined in (I6.14p . which are zero in the light of Lemma [3j 

Proof for ( 16.10dl) : Using (16.1 Pel) to eliminate the term (a — b)fg from ( I6.10dl) . we get 

J[f - (a + 5)g\ - J[f- (6 + 6)g] + (a - 6)7^ = 0. (6.18) 

We fix k = 3 and c = S in (jCTTjl . and use the formulas fjOldL (ICHeT) . dCH and fjClhl) . 
Then it turns out that 

7[/ - (a + 5)?] - Jif - (6 + 5)g] + (a - b)Jg 
= (a + 5) 2 Z X - (6 + 5) 2 Z 2 - (a + 5) 2 (6 + 5) 2 Z 3 

with Zj defined in (I6.15p . which are zero in the light of Lemma [3j 
Thus we have completed the proof for all bilinear equations. 

□ 
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7 Ql with power background 

7.1 Background solution with T(x) = —x + c 

With fixed point denned by T(x) = —x + c we use the reparameterization 

p = i r (l - cosh(a')) = -\r{l - a) 2 /a, (7.1) 
q = i r (i _ cos h(/5')) = -|r(l - /3) 2 //?, (7-2) 
x nm = Ae w + Be Vnm + c/2, where = 5V/16, (7.3) 

and this leads to equations that factorize as in (15. 4p . Thus we get power type background 
solutions for u: 

u = \c + Aa n f3 m + Ba~ n p- m , AB = 5V/16. (7.4) 

Here c is related to translation freedom and r to scaling freedom; in the following we take 
c = 0. 

7.2 1-soliton solution for Ql s 

In order to derive the 1SS we use the BT (16.41) with the seed solution (17.41) with 

u — Uq + v, (7.5) 
where u is the bar-shifted background solution (17.41) : 

u = Aa n p m K + Ba- n p- m K-\ AB = <fV/16. (7.6) 
and k is defined through 

x=-^(l- K ) 2 /*, (7.7) 

If we now define 

U n , m = Aa n (3 m - Ba- n (3~ m (7.8) 
then it is straightforward to derive (12.111) from (16 .4p with 

r(l- a){\- k){\- aK) _ r(l - q)(l - /c)(q - k) 

m r(l-fl(l -«)(!-/?«) n _ r(l-/?)(l-«)G9-«) 

On the basis of this we define p as usual by 

(S\ n (T\ m (l-otK\ n (l-(3K\ m 

p ^=[a) UJ Po ' 0= U^J b^J Po ' ' (7 - u) 

where po,o is some constant. Then it follows that 

!~« 2 77 n 
i" Pn,m 
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and finally we obtain the 1-soliton for Ql 5 : 

u l S i = UQ + V n:m (7.13) 
Aa n f3 m (K + «-Vn,m) + Ba- n f3- m (K~ l + K Pn , m ) 



1 ~\~ pn,m 

A'a n (3 m {l + K- 2 Pn , m ) + B'a- n f3- m {\ + /t 2 pn, m ) 



(7.14) 
(7.15) 



1 ~\~ Pn,m 

where A'B' = AB = 5 2 r 2 /lQ. 
7.3 Bilinearization 

In order to get p of (17.1 ip into a nicer form we use the Mobius transformations 

a — c b — c k — c rc 2 rc 2 

a + c o + c fc + c cr — c z cr — cr 

which leads to the canonical form 
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. a + A;\ n (b + k\ 
Pn,m = r 7 — r Po,o- (7.17) 



a — k J \b — k 

With the above Mobius transformations Ql 5 (16.11) can be written as 
Q l< = _ c 2)( u _ _ g) _ ( a a _ c 2)( u _ 5)( S _ ~) _ ^g^lj^ = o, (7.18) 

Using (17.17P and rearranging the parameters A and B we write the 1SS (17. 15ft as 

^n,m,l) +Ba P rf,(n,m,l) ' — ^ r /16, (7.19) 

where ^ is as in (12.201) . 

This structure motivates us to bilinearize Ql 5 through the following transformation, 



7 / 

u™? = Aa n (3 mJ j + Ba- n (3~ rn J, AB = 5V/16, (7.20) 



where / = \N — 1L , (z/ = 1,2, or 3) with entries (I2.20p and the bar-shift is the shift in 
the third index I, as discussed in Sec. 12.51 

The solution (I7.20p with r = 1 is the same as (1-5.11) with B — 0. In fact, by a 
comparison ( 17. 16ft with (1-3.17) one has 

p(a) = a~ n p- m , p(-a) = a n [3 m ; 



28 



and substituting pi in (1-2.2) by (15.161) we find from (1-3.17) that 

7 / 

(1 + 2aS(-a, -a)) x - a) 2 ) = L, (1 - 2aS(a, a)) / (U(kj - a) 2 ) = j. 

Surprisingly Ql* 5 can also be bilinearized with the same bilinear equations as H3 5 , 
namely (I5.17P , and we have already shown in Proposition S] that / = | JV — 1 1 solves the 
Bt in fl5TTj) . 

The representation of Ql 5 (17.181) in terms of Bi is as follows: 

rn« - a 4 " +2 /? 4m+2 ( a + c) 2 (b + c) 2 A 2 P 1 - a 2n (3 2m 5 2 /4P 2 - (a + c) 2 (& + c) 2 B 2 P l 

a 2n + x (3 2m+x {a + c) 2 (6 + cfffff 

YY-XX, X = B 1 -2cff, Y = B 2 -2cff, 
_( a + c )( a _ c )( 6 + c )2 p _ 4c 2 77/7) 

_( a + c )( a _ c )( 6 _ C )2(xf - 4c 2 / J??) 

+4c 2 (6 + c)(6 - c) (XX - Ac 2 fffJ) 

+( 6 + c )( 6 _ c )( a + c )2^yf - 4c 2 777/ ) 

+ ( 6 + C )( 6 _ c )( a _ c )2 (yF - Ac 2 fjjj) 

-Ac 2 (a + c)(a - c) (f? - 4c 2 ////) . 

8 Conclusions 

In this paper, companion to [1] we have analyzed the soliton solutions to the models 
HI, H2, H3, and Ql in the ABS list |2J of partial difference equations. Our method 
is constructive, progressing in each case from background solution to 1SS to NSS and 
bilinearization. 

Our approach is fairly algorithmic, and as such we hope that it will be usable also for 
other models with multidimensional consistency. One interesting feature is that in some 
cases we need several bilinear equations, some of which seem to have the same continuum 
limit. This is a remainder that there are several ways to discretize a derivative. 



where 

Pi 

P2 
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A Casoratians and formulas 

We define Casoratians 

/ M = PV^1| H , g M = \N^2,N\ M , h M = \N^2,N + l 
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where the matrix entries are given by the function 

^(n, m, I) = gt(c + h)\a + h) n (b + k^ + g-(c - h)\a - h) n (b - k t ) m 
We introduce notations 

r V)i = o£ - A;?, r i/ = Diag(r i/) i,r i/i 2, ••• ,T VjN ), (1/ = 1,2,3), 

where 

ol\ = a, CH2 = b, as = c, 

o 

and define the operator E u as 

E^ = T- 1 E^, (i/ = l,2,3). 
The column vector if) satisfies 
(a, - a K )iP = (£" - E K )ip, 
a v )E»E v i) = (E» + E v )ip, 

E^ = [E K - (a„ - a K )(E K ) 2 + (t*„ - a K ) 2 ^(E K ) 2 ]^, 

= - K~ a K )(£ K ) 2 + (a„ - a K ) 2 (£ K ) 3 - (a„ - a K ) 3 E^E K ) 3 ]^A.5d) 

For the Casoratian f, . , we have 

^ [k] ' 



(A.l) 

(A.2) 
(A.3) 

(A.4) 

(A.5a) 

(A.5b) 
(A.5c) 



- K - a K ) N - 2 E,f [K 


= \N-2,E^(N-2)\ [k] , 




(A.6a) 


K - a^EJ^ 


= \N^2,E^{N-1)\ [kV 




(A.6b) 


K + a K ) N ~ 2 E»f [K 


= \T,\\N~^2,E^(N-2)\ 




(A.6c) 


K + a K ) N ~ l E»f [K 


= \T,\\N^~2,E^(N-1)\ 




(A.6d) 




= \Kv 




(A.6e) 




= |-l,jV^2| w , 




(A.6f) 


K - a K ) N ~ 2 E^f [K] 


= \N-1,E^(N-1)\ [k] , 




(A.6g) 


K - a K ) N ^E^EJ [K] 


= \-l,N^3,E^(N-2)\ [ 




(A.6h) 


(a, + a K ) N ~ 2 E»E K f [K] 


= \t,\\n^i,e^(n-i)\ m 


, (/i ^ re), 


(A.6i) 


(a, + a K ) N - l E»E K f [R] 


= |r M || -l,N^3,E»rl>(N- 


2)I M , 0^«); 


(A.6j) 



(a M - a*) (a,, - a K ) {a v - a K ) E^E u f [K] 

=\N^3, E^(N - 2), E^(N - 2)| M , 
K + «„)(«„ + a K ) N ~ 2 {a v - a K ) N ~ 2 E^E v f 



(A.6k) 



|rV||iV - 3, E v rl>{N - 2), E^(N - 2)1,, (/x ^ 



(A.61) 
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= \N - 2, EMN - 1), E^(N - 1)| M , (A.6m) 
+ c^)K + a K ) N ~\a u - a^) N ~ 2 E»E u E K f [K] 

=\r^\\N^2,E u ilj(N - 1),E^(N - 1)\ [k] , (ft + «). 

(A.6n) 

For o. . we have 

- (a ll -a K f- 2 E ll [g [fi] + (a^-a^fj = \N^3, N - 1, E^{N - 2)| M , (A.7a) 
K+a K ) 7V - 2 ^[^ ] -K+a K )/ H ] = |r M ||iV^3,iV-l,^(iV-2)| M , ^ «XA.7b) 

For /i w , 

-(a M -aX"X^ M +K- "«)» M ] = |A 7r 3,iV,^(iV-2)| H , (A.8a) 
(a, + a/- 2 %,-(«/. + ««)ff M ] = MA^, N, E^(N - 2)| [k] , «)(A.8b) 
These formulas, ( 1A.5I) - (1A.8I) . are valid for fj,,u,K— 1, 2, 3 except when otherwise indicated. 

B Proof for the formulas in Appendix [A] 

The formulas ( 1A.5I) can directly be verified. In ( 1A.6I) - (1A.8I) . those only containing down 
shifts Efj, and E v can be proved by using the relations (1A.5I) . Since a, b and c appear in 
ipi equivalently, we take fi = 1, v = 2 and k = 3 in the following proof. 
From (IA.5al) we have 

(a-c)f(l) = ^(0-^ + 1), (B.l) 
(b-c)i>{l) = + (B.2) 

(0-6)^(0 = ^(i)-i(O- (B.3) 



Using fIB.ll) for / . we first find 



(3) 

(a - c)f [3] = | (a - c#(0), #(1), . . . ^(N - 1)| [3] = |^(0), #(1), . . . ,^(A - 1)| [3] 
and repeating this N — 1 times we get 

(a - c^Vp, = |V(0), • • • ,^(A - 2),^(JV - l)| m = |A^2, ^(A - 1)| 



This is formula ( 1A.6bl) . Again using ( IB. 1}) and expressing ip(N — 1) through ^(A — 2) 



and ^(A — 2) in the above formula yields 



(a - c)"- 2 / m = | A - 2,f(N - 2)| [3]) (B.4) 
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which is (lA.6a|) . For / we first using (1B.2j) from the above formula to obtain 



(6 _ c )^( a _ c) ^ f[3] = [iV 3, ^(iV - 2), V(iV - 2) | 



then from (IB .3|) we get 



(a _ 6 )( a _ c )^(& _ ^a-^ = |iv _ 3 , V»(JV - 2),^(JV - 2)| 



[3]' 



(B.5) 



which is (IA.6kp . 

For g = |JV^2,JV| [3] we find first 



,,, - N-3.v(N-2),v(N)\ 3 .. 
Using (IA.5cl) to re-write the last column we get 

(a - c) N - 2 g [3] = -\N^3, N - l,f(N - 2)| [3] + (a - c)\N^2,f(N - 2)| 
which can also be stated as 

- (a - c) N - 2 [g VA] + (a- c)fj = \N^3,N- 1^(N - 2)| 



[3]' 



[3] ' 



(B.6) 



(B.7) 



i.e., ( IA.7al) . Formula (1A.8aj) can be derived similarly by using (1A.5dj) . 

To show how to derive those formulas with up shifts in (1A.6l) -( TA~8l) . we introduce two 
auxiliary functions: 

II : <f>i(n,m,l) = Qj{c + k i )\a + k i ) n {b-k i )- m + QT(c-k i )\a-k i ) r \b + k i )- m , (B.8a) 
III : ufa m, I) = et(c + h) l (a - k t )~ n (b + h) m + 6 ;(c - h) l (a + h)- n (b - h) m . (B.8b) 



Casoratians /, g and h (w.r.t. bar-shift) with column vectors 



) T and 



UJ 



(wi,--- ,u N ) T are denoted by / II[3] , g n[3] , h n[3] and /, 



They are related to / [3] , g [3] and h [3] through 
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III[3]' i>III[3]' III[3]' 



f 13] I "^2 | / n [ 3 ] 

|r 2 l rri ^ II[3] 



9m 



1-^1 1 /lll[3]' 
l-^ll ^III[3]' 



[3] 



\r 2 rh m = ir^h 



III[3] ' 



respectively. 

(B.9a) 
(B.9b) 
(B.9c) 



which follow from ip = r™0 = T™uj. 
6 satisfies 



(a — c) (f)(1) 
(b + c)$(l) 
(a + b) 



0(0-^ + 1), 

0(0 + 0(0- 



(B.10) 

(B.ll) 
(B.12) 
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Similar to the previous case of ip, using these relations we can get 



(b + cf- 


~ 2 f 

Jll[3] 


= \N- 


~2,0(iV- 


2) II[ 3 ] ' 


(B.13a) 


(b + c) N ~ 


Jn[3] 


= |JV- 


~2,0(iV- 


■'") II[3] ' 


(B.13b) 


a + b)(a-c) N ' 2 (b + c) N - 


in[3] 


= 




2),0(iV-2)| II[:j] ; 


(B.13c) 



(b + c) N - 2 J im = \N — 1, <p(N — 1)| II[3] , (B.13d) 
_ (a + 6 )( a _ c )iv-2 (6 + c) iv-2j ii[3] = |Arr^,0(iV-l),0(iV-l)| n[3p ( B .13e) 

(6 + c)^% 3] = | -1,JV^3, J(JV-2)| II[3] ; (B.13f) 



(& + c )^[ ?ii[3] _( fo + c ) /ii[3] ] = |iV-3,iV-l,0(iV-2)| n[3] , (B.13g) 
(b + c) N -% m -(b + c)g im ] = \N^3,N^(N-2)\ U[3] . (B.13h) 

Now, noting that / [3] = |r 2 | ?n / II[ 3 ] and ip — r™0, we have 
(b + c) N -% = (b + c) N ~ 2 \T 2 r +1 f m] = \T 2 \\m,---,m-2),E^(N-2)\ [3] , (B.14) 

where the operator E 2 is defined as ( 1A.4I) . The above formula is just (1A.6cl) for fi = 2. 
From the rest of (IB. 131) we can get other formulas with an up-hat shift in ( 1A.6I) - (1A.8I) . 
where for (1A.61I) and ( lA.6nl) (p = 2,v = 1) we need to use (lA.5bl) . Besides, if we take a 
down-hat shift for (IF3. 14[) and rewrite the last column by using flA.5b|) (with fi — v — 2), 



then we get 

26(6 + c) N ~ 2 (b - cf- 2 f [3] = \T 2 \\N^3, f(N - 2), E\{N - 2)| M , 

i.e., the formula flA.611) of the case fi = v = 2. 

Next, using u we can derive these formulas with an up-tilde shift in (lA.6j) - f)A.8j) . uj 
satisfies 

(a + c)w(l) = u(l)+u{l + l), (B.15) 

(b-c)uj(l) = - uj(1 + 1), (B.16) 

(a + b)u(l) = w(l)+u(l), (B.17) 

and from these relations we have 



\N-2i 



N-2 7 

J III[3] 


= |jv^ 


2, 


E(iV- 


2) m[3]' 


(B.18a) 


N-1 7 

J III[3] 


= 1^ 


2, 


w(JV- 


1) rap] ' 


(B.18b) 


Af-2 7 

I HI[3] 


= |jv^ 


3. 




2),uJ(JV-2)| III[S] ; 


(B.18c) 
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(a + c) N - 2 J um = \N - - l)\ nm , (B.18d) 

(a + b)(b-c) N - 2 (a + cf- 2 J um = \N^2,uj(N-l)MN-l)\ um , (B.18e) 

(a + cf- 1 /^ = \-l,N^~3MN-2)\ m[3] ; (B.18f) 

(a + c) N - 2 [g lim -(a + c)f lim } = \N^3,N - 1,Z(N - 2)| m[3] , (B.18g) 

(a + c) iV - 2 [V [ 3 ] -(a + c)^ III[ 3 ] ] = 1^^3^,2(^-2)1^. (B.18h) 

o _ _ 

By using operator E , these formulas will generate those of (1A.6fl - flA.8h with an up-tilde 
shift. 

Thus we can get all formulas for /, g and h given by (1A.6h - f1A.8jl . 



C Formulas for f = \N - 1| [3]J g=\-l,N- 1| 

For the Casoratian f. , , we have 

= E3/-K-«3)|^(-l),iV^2| [31) (Cla) 

E^E 3 / = f-(a,-a 3 )g+(a,-a 3 ) 2 \E^(-l),N^l\ m , (C.lb) 

£^ = |^(-l),iV^| [3] -K-«3)|^(-l),-l,iv^2| [ 3 ] , (C.lc) 

-E»f = E 3 f - (a, + a 3 )\E^(-l),N^2\ m , (C.ld) 



:-i) N 



E^E A f = f + ( afl + a 3 )g-(a, + a 3 ) 2 \E^(-l),N-ll 3] , (C.le) 



A' I 



= -l^(-i) ) ^^l [ a ] -K + «3)l^(-i),-i,A rr 2| [3] ;(c.if) 

(a - 6)/ 

= ( fl _ b )f _ ( a _ c) V(-l), iV^2| [3] + (b - c) 2 \f(-l), N^2\ [3] 
+ (a - C ) 2 (6 - c)|^(-l), -1, N^2\ [3] - (a - c)(6 - c) V(-l), -1, iV^2| [3] 
+ ( a _ C )2( 6 _ c )2|^(-i) s ^(_i) s aT^2| [3] , (C.lg) 

a ~ b J 



|ri||r 



2 



( a _ 6 )/ _ ( a + C ) 2 |EV(-1), iV - 2| [3] + (6 + c) 2 |EV(-l), N - 2| [3] 

_ ( a + c )2( 6 + C )|^V(-1), -1, A^2| [3] + (a + c)(6 + c) 2 |E 2 ^(-l), -1, aT^2| [3] 

_ ( a + C )2( 6 + C ) 2 |E 2 ^(-1), EV(-l), N^2\ w (C.lh) 
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D Proof for the formulas in Appendix [C 

Noting that 

$(l + l) = 1>(l)-(a-c)f(l), (D.l) 

we have 

/ =1^(0)^(1), ■■■,f(N- 3), f(N - 2),f(N - 1)| [8] 
=\f(0),f(l), ■■■,iP(N- 3), f(N - 2), ^(N - 2) - (a - c)f(N - 2)| [3] 
= |^(0), ^(1), • • • , 4>(N - 3), ^(iV - 2), ^(iV - 2)| [3] 
= |^(0), • • • , V>(iV - 3), ^(iV - 3), V(iV - 2)| [3] 



= |^(0), ^(0), • • • , V(iV - 4), ^(iV - 3), V(iV - 2)| [3] 
=|-0(— 1) - (a - c)^(-l), V(0), • • - , V(iV - 4), ^(iV - 3), V(N - 2)| [3] 
=/-(o-c)|^(-l),iV^2| M . 
This is (IC.lal) for jj, — 1. In a similar way we can prove (IC.lbl) . (IC.lcl) and (C.lg) 



To prove (IC.ldl) . ( IC.lel) . (IC.lfl) and (IC.lhj) . we consider the Casoratian type determi- 
nants / = \N — 1| [3] , g — | — 1, N — 1| [3] with the following entries 

IV:<7i(7i,m,0 = f ? +(c + ^)'(a-A; l )- n (6-A; 4 )- m + ^( c -A; l )'(a + A; J )- n (6 + A; i )- m , (D.2) 
and denote such / and g by / IV[3] and g IV[3] . Noting that a satisfies 

(a + c)a(l) = a(l) + a(l + l), (D.3) 

(b + c)a(l) = a{l)+a{l + l), (D.4) 

(a-b)a(l) = a{l)-a(l), (D.5) 



we can have 



-l) N E»f lvm = E 3 f m3] -(a, + a 3 )\E»a(-l),N-2\ msv (D.6a) 



'IV[3] ^JJTVIS] 1 "Jyi-^ ^ V V "llV[3]' 

(-l) N E»E*f IV[3] = / IV[3] + (a„ + a 3 )<? IV[3] - K + a 3 ) 2 \E^(-l), N - 1| IV[3] , (D.6b) 

-K + as)|SV(-l), -1, iV^2| IV[3] ; (D.6c) 

( a " 6 )7 IV[3] 

= (a - 6)/ IV[3] - (a + c) 2 |?(-l), jV-2U m + (6 + c) 2 |a(-l), A^2| IV[3] 

_ ( a + c )2 (6 + c)ISr(-l), -1, iv^2| IV[3] + (a + c)(b + c) 2 |a(-l), -1, N^2\ W[3] 

- (a + c) 2 (6 + c) 2 |a(-l), 3F(-1), iV^2| IV[3] . (D.6d) 

Then, using the relationship / I[3] = iT^T^f g m = {^{T^g^ and ^ = T^Tfa, 
we can finally derive the formulas (1C. ldf) . (IC.lejl . (IC.lfl) and (1C. lh|) . 
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